function [ ] = plotEnergy( tri, rx, ry, rho )
%PLOTENERGY Summary of this function goes here
%   Detailed explanation goes here

% cons(:, 1) = domain2D.e1;
% cons(:, 2) = domain2D.e2;
% 
% dim = size(rx, 1); 
% 
% px = [domain2D.vx; rx];
% py = [domain2D.vy; ry];

px = [rx];
py = [ry];

n = size(px, 1);

c = zeros(n, 1);

c = rho.^3;

figure;

% dt = DelaunayTri(px,py,cons);
% io = dt.inOutStatus();
% tri = dt(io, :);

trisurf(tri, px, py, c, c, 'facecolor', 'interp', 'edgecolor', 'none');

colormap(jet(256))

view(2);

axis ( [ -1, 1.0, -1, 1.0] );
caxis( [ 0, 220000] );
axis square;

colorbar;


end

